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Abstract. We present broad band power spectra of variations of SS433 in radio, optical and X-ray spectral bands. 
We show that at frequencies lower than 10~^ Hz the source demonstrates the same variability pattern in all these 
bands. The broad band power spectrum can be fitted by one power law down to frequencies ~ 10~^ Hz with 
flattening afterwards. Such a flattening means that on time scales longer than ~ 10^ sec the source variability 
becomes uncorrelated. This naturally leads to the appearance of quasi-poissonian flares in the source light curve, 
which have been regularly observed in radio and optical spectral bands. The radio flux power spectrum appears to 
have a second break at Fourier frequencies ~ 10~^ Hz which can be caused by the smearing of the intrinsic radio 
variability on timescale of the light-crossing time of the radio emitting region. We flnd a correlation of the radio 
and optical fluxes of SS433 and the radio flux is delayed by about ~ 2 days with respect to the optical one. Power 
spectra of optical and X-ray variabilities continue with the same power law from 10~^ Hz up to ~ 0.01 — 0.05 Hz. 
The broad band power spectrum of SS433 can be interpreted in terms of self-similar accretion rate modulations in 
the accretion disk proposed by Lyubarskii (1997) and elaborated by Churazov et al. (2001). We discuss a viscous 
time-scale in the accretion disk of SS433 in implication to the observed broad band power spectrum. 
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1. Introduction 

SS433 is a high-luminosity massive X-ray binary sys- 
tem with steadily precessing accretion disk and jets 
( |Margon 1984[ see IFabrika 2004 I for a recent review). 
Most of the system luminosity (~ lO''" erg/s) is emit- 
ted in optical and UV spectral bands (Cherepashchuk 
et al., 1982; Dolan et al., 1997). The compact source in 
this binary system, probably a black hole, accretes the 
material from the companion late A-supergiant (Gies et 
al. 2002; Hillwig et al. 2004; Cherepashchuk ct al. 2003; 
2004) at a highly super-Eddington rate. Mildly rclativistic 
{v 0.26c) baryonic jets are launched from the innermost 
parts of the supercritical accretion disk. 

The system demonstrates several periodicities 
(Cherepashchuk 2002): the precessional one (~162 days), 
the orbital one (~13 days), and the nutational one 6 
days). Studies of these periodicities made it possible 
to tightly constrain some of the system's parameters 
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(Eikenberry et al. 2001; CoUins & Scher 2002) including 
the binary system inclination angle {i ~ 78°) and the jet 
precession angle 20°). 

X-ray emission produced in the innermost regions of 
the accretion flow is completely screened from the line of 
sight by the optically and geometrically thick inner disk 
and outflowing wind. The observed X-ray flux from SS433 
can be modeled by thermal bremsstrahlung emission from 
a ~ 20 keV plasma moving in the jets (Watson et al. 1986; 
Kotani et al. 1996). A detailed analysis of high resolution 
X-ray spectra taken by CHANDRA (Marshall et al. 2002) 
basically confirms previous findings. The recent analysis 
of hard X-ray spectra of SS433 obtained by INTEGRAL, 
however, indicates that a mixture of plasmas with different 
temperatures and geometries could be present in the jet 
region (Cherepashchuk et al. 2003; 2004). 

The source demonstrates variability in all spec- 
tral bands and on various time scales (see e.g. 
IFiedler et al. 19871 Goranskii et al. 1987, Zwitter et al. 
1991, Kotani et al. 2002). However, thus far no systematic 
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studies of aperiodic variability of SS433 at different wave- 
lengths have been carried out. Previous works were mainly 
focused on the analysis of periodic variabilities found in 
the source. SS433 has been regularly monitored at differ- 
ent wavelengths over the last decades and appears to be 
also well suitable for analysis of aperiodic variations. The 
long-term optical observations indicate that the source 
demonstrates on average a very stable behavior with no 
drastic changes in its activity pattern (Eikcnberry ct al. 
2001; Goranskii et al. 1998). It is the only known per- 
sistent galactic microquasar with nearly constant kinetic 
power in jets (Fabrika 2004). 

An enhanced radio flux is observed from the source at 
some periods. These were classified as active states on top 
of the quiescent (passive) states (Bonsignori-Facondi et al. 
1986; Fiedler et al. 1987; Trushkin et al. 2003). However, 
on the long range, it can be shown that the appearance of 
the active radio states in SS433 is aperiodic process and 
can be well described in the language of power spectra. 

In this paper we for the first time attempt to system- 
atically study the aperiodic variability of SS433 by con- 
structing power spectra of its flux variations in radio, op- 
tical and X-ray spectral bands. We found that variabilities 
at different wavelengths can be described by one power law 
spanning several orders of magnitude in frequency, which 
is suggestive of their common nature. Wc discuss the pos- 
sible origin of the obtained power spectra as being due 
to the mass accretion rate modulation in the underlying 
accretion disk. 



2. Observations and data analysis 

To estimate power spectra of the source variability in dif- 
ferent spectral bands we have combined all data avail- 
able to us. As the data were obtained by different instru- 
ments, are not evenly spaced and have significant time 
gaps, we can not use the simple discrete Fourier trans- 
form, which is usually applicable for evenly spaced data 
(e.g [Leahy et al. 1983| ). 

The optical and radio data are analyzed by means 
of periodograms that provide estimations of the Fourier 
amplitude at a given Fourier frequency / in the form 
dPeeming 1975llLomb 19761 |Scargle 1982| ): 
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Here f{ti) is the flux measurement at a given time U with 
subtracted average flux value, N is the number of mea- 
surements in the set. To obtain the fractional rms squared 
normalization we calculate the power P{f) in the form: 



Pif) 



2T 



HI)\ 



where T is the time span of the set {T ~ t^q ~ ti). The 
power spectrum calculated with this normalization has 



one important property: The integral of the power spec- 
trum from /i to /2 gives the square of fractional rms vari- 
ability which is present in the lightcurve of the source on 
time scales from /j"^ to f^"^. 

The longest time series were divided into a few parts 
and the power spectrum was evaluated for each of them. 
Estimates of the power of different parts lying at close 
Fourier frequencies were averaged and the standard devi- 
ations of these values were calculated. 

For all our data sets we have made simulations in order 
to check the effect of uneven sampling. Simulations show 
that this effect distorts the power spectrum in the form of 
a power law [P oc /~") at the low frequency end. but not 
very significantly. In most cases such a distortion does not 
exceed the uncertainty in the power estimates. 



2.1. SAI database of observations of SS433 

The optical broad band power spectrum of SS433 was 
calculated using a data-base of optical V-band photom- 
etry of SS433 obtained in 1979-1996. The data-base con- 
sists of 2200 individual nights of observations collected at 
Sternberg Astronomical Institute (jCoranskii et al. 1998|l . 
In this data-base all published V-band observations are 
reduced to one photometrical system. The average uncer- 
tainty of the photometrical data in individual nights is 
5-10%. The data-base also includes series of photomet- 
ric observations taken during several consecutive nights 
(for example, Goranskii et al. 1987; 1997), with each night 
containing from several to a few tens individual observa- 
tions. The contribution of measurements uncertainties to 
the obtained power spectrum (the white noise component) 
have been subtracted. Besides, we have removed from the 
power spectrum the peaks due to coherent modulations 
with orbital, nutational and processional periods. 



2.2. RTT150 

The high-frequency power spectrum of optical fluctuation 
was calculated using observations of SS433 performed with 
1.5-m Russian-Turkish Telescope (RTT150) at TUBITAK 
National Observatory (TUG), Bakyrly mountain. The 
observations were carried out in September 2004 dur- 
ing seven nights under clear sky but poor seeing (« 
2") conditions. In our analysis, we used points out of 
the orbital eclipse. A low readout noise back-illuminated 
2x2K Andor Technologies DW436 CCD mounted in 
F/7.7 Cassegrain focus of the telescope was used. The pho- 
tometric R filter was utilized because of lower Galactic 
absorption in this filter and larger magnitude of SS433 
(R« 12) than in blue bands. The method of fast photom- 
etry we used is described in some detail in Revnivtsev et 
al. (2004) 
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2.3. X-ray data: RXTE/ASM and EXOSAT/ME 

In order to construct the broad band power spectrum 
of variability of SS433 in X-ray band we have used 
two sets of data. Until now, the longest time-series of 
SS433 in X-rays is provided by the All Sky Monitor 
(http://xte.mit.edu/ASM_lc.html) onboard the Rossi X- 
ray Timing Explorer IjBradt et al. 1993)l . 

The power spectrum calculated from the observed X- 
ray ASM light curve should have the statistical noise 
contribution (a constant component in the power spec- 
trum of the source). The value of this constant compo- 
nent in the power spectrum in principle can be calcu- 
lated theoretically by using the values of statistical un- 
certainties quoted in the light curve. However, the er- 
rors quoted in the ASM light curves often underesti- 
mate the real uncertainties of the measured fluxes (see 
e.g. [Grimm, Gilfanov, fc Sunyaev, 2002| ). Therefore in or- 
der to estimate the contribution of the noise component we 
have studied the power spectra of persistent weak sources, 
the supernova remnants Gas A and Puppis A. The ob- 
tained noise component after the appropriate renormal- 
ization was then subtracted from the power spectrum of 
SS433. 

To construct the X-ray power spectrum at frequencies 
> 2 X 10"^ Hz we have used data of EXOSAT/ME obser- 
vations. The EXOSAT observatory had orbit much longer 
than that of RXTE thereby providing much longer unin- 
terrupted sets of observations. Preprocessed lightcurves 
of SS433 were taken from the HEASARG archive at 
GSFG (http://legacy.gsfc.nasa.gov). Only observations 
performed at precessional phases 0.8 < tp < 1.2 were con- 
sidered. At these precessional phases the X-ray flux from 
jets of SS433 is maximal. Power spectra were constructed 
using the standard task of FTOOLS 5.2 package. 

2.4. Radio observations 

The longest available set of radio flux from SS433 
at frequencies 2.25 GHz and 8.3 GHz is provided 
by the GBI monitoring (http://www.gb.nrao.edu/ fg- 
docs/gbi/gbint.html). We have used measurements at 
the frequency 2.25 GHz because these have higher sig- 
nal to noise ratio than those at 8.3 GHz. The contribu- 
tion of statistical noise was subtracted from the obtained 
power spectrum. The available data allowed us to esti- 
mate the power at frequencies from ~ 5 x 10~^ Hz to 
~ 2 X 10"^ Hz and at 10~^ Hz. Before the binning 
of the power spectrum we have removed the peak of co- 
herent variations with the nutational period (~6 days, 
[Trushkin, Bursov, fc Smirnova "20011 ). 

In order to extend the power spectrum of radio vari- 
ability of SS433 to higher frequencies we have included 
the estimate of the SS433 variability on time scales of 
tens and hundreds of seconds made by Rand et al. (1988). 
Note that these measurements were done at the frequency 
0.4GHz (different from 2.25 GHz for the GBI measure- 
ments), so the relative normalization of the power esti- 



mates might be different. However, comparison of avail- 
able long time scale measurements of radio flux of SS433 
at 0.4 GHz (e.g. |Bonsi gnori-Facondi et al. 1986| with that 
at 2.25 GHz tells us that the normalization correction fac- 
tor between variabilities at these two frequencies does not 
exceed ~ 2 — 3. 
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Fig. 1. The power spectrum of the SS433 variability in 
different spectral bands. The optical power spectrum is 
shown by crosses, the radio power spectrum is shown by 
open circles, the X-ray power spectrum is represented by 
filled circles. No special rescaling is done. The dashed line 
shows the fit in the form P(/) oc (1.0 + (///broak)")"\ 
where a = 1.5 



3. Results 

The resulting broad band power spectra of aperiodic vari- 
ability of SS433 in different spectral bands is shown in 
Fig.n Two most important features of the power spectra 
are clearly seen. 

— First: stochastic flux variabilities in all three energy 
ranges (radio, optical and X-ray) follow very similar 
patterns. The power spectra in these ranges are flat- 
tened in Fourier frequencies ^ 10~^ Hz. 

— Second: the power-law part of the optical and X-ray 
power spectra (which can be traced up to Fourier fre- 
quencies of the order of 0.01-0.1 Hz) spans over more 
than five orders of magnitude in frequency. The power 
spectrum in radio is likely to have a break between 
IQ-^ Hz and lO"'' Hz. 
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3.1. Self-similar power spectrum - is it produced by 
accretion disk? 

The self-similar behavior of variability of sources is ob- 
served not for the first time. Decades ago it was noted 
that at low Fourier frequencies (< 0.01 — O.lHz) X- 
ray binaries demonstrate a self-similar power law compo- 
nent in their power spectra (e.g. Ivan der Klis et al. 19871 
[Hasinger fc van der Klis 1989| ). A detailed study of vari- 
ability of the black hole binary system Cyg X-1 in the 
soft spectral state (when the optically thick accretion 
disk presumably extends down to the last stable orbit 
around the black hole) have demonstrated that this power- 
law component can extend over four orders of magnitude 
in Fourier frequency (jChurazov et al. 2001(1 . An extensive 
study of low mass X-ray binary systems revealed that such 
power-law behavior is observed in practically all LMXBs 
(Gilfanov & Arefiev 2005). 

It was suggested that the observed self-similar variabil- 
ity in X-ray binary systems is a result of mass accretion 
rate variations in the accretion disk around the compact 
object ( |Lyubarskii 1997[IChurazov et al. 2001L Gilfanov & 
Arefiev 2005). The observed variability in X-ray binaries 
at all Fourier frequencies is rather strong (a significant 
part of the total luminosity varies), which implies that in- 
ner regions of the accretion disk, where the most of the 
energy release takes place, are involved. These regions, 
in turn, respond to the accretion rate modulations which 
have diffused there from the outer parts of the disk. If 
modulations of the accretion rate occur at any given ra- 
dius of the accretion disk, imprinting the characteristic 
frequency to the inflowing matter, one can readily get 
the self similar power spectrum of the form P oc /~", 
where a 1 — 2 for a wide range of accretion disk models 
ULyubarskii 1997| ). 

3.2. Viscous time scale in the accretion disk 

Such a power-law shape of the accretion disk power 
spectrum is expected to extend up to time scales on 
which the disk variability becomes uncorrelated, i.e. 
up to the longest time scale of the accretion disk 
IjGhurazov et al. 20011 Gilfanov & Arefiev 2005). From 
theoretical point of view this time scale is the time 
of viscous diffusion of matter from the outermost re- 
gions of the accretion disk to the innermost parts. 
For example, in the framework of standard optically 
thick accretion disk fShakur a fc Sunyaev 1973| see also 
|Lynden-Bell fc Pringle 1974; ): 

(H\-^ Tk{R) 

where H is the hydrostatic height of the disk at radius i?, 
Tk{R) is the Keplerian rotation period at this radius, a 
is the viscosity parameter. 

At frequencies below t~-^^^ the variability becomes 
uncorrelated and consequently the power spectrum should 



flatten. As mentioned above, this indeed was found to be 
the case for all LMXBs (Gilfanov fc Arefiev 2005). The 
uncorrelated variability at frequencies / < t~-^^ should ap- 
pear on the light curves of the source in different energy 
bands as a set of random flares with durations around ~ 
tvisc- Such "flares" have been actually often observed in ra- 
dio and optical spectral bands, which led to the classifica- 
tion of the so-called "active" and "passive" states of SS433 
(e.g. [Bonsignori-Facondi et al. 1986| llrsmambetova 19971 
IFabrika^0^^r ~Such a flaring behavior in reality ap- 
pears to be a natural observational appearance of the 
uncorrelated variability on time scales larger than the 
longest (viscous) time scale of the accretion disk. Note 
that a very similar phenomenon (but on much shorter 
time scales) has been observed many times from the best 
studied black hole candidate Cyg X-1, whose variability 
can be described as a series of random shots (the "shot 
noise" model, e.g. ITerreU 19721 IVikhlinin et al. 19941 
[Negoro, Miyamoto, fc Kitamoto 1994| ) 

The break in the power spectrum of SS433 at / ^ 
1 — 2 X 10^^ Hz corresponds to a time scale of 50-100 
days. The ratio of this time scale to the orbital period 
13.6 days in SS433 is ~ 5 — 10, which is in a reasonable 
agreement with results obtained by Gilfanov & Arefiev 
(2005) for LMXBs. ^ Such a low value of the ratio may 
indicate that the outer accretion disk is relatively thick 
H/R 0.1 — 0.2 (see also the discussion in Gilfanov fc 
Arefiev 2005). In the case of SS433 the thickness of the 
outer parts of the disk should be at least not smaller than 
that found in LMXB systems because of the very high 
accretion rate. So in principle for a low alpha-parameter 
~ 0.1 the viscous time scale of the accretion disk in SS433 
evaluated from the above formula can be comparable to 
the value derived from the power spectrum analysis. 

Note here that the ratio of orbital to the break fre- 
quency of the order of 5-10 may indicate that the mass 
ratio in the binary system SS433 is larger than 0.35 
(see results of Gilfanov fc Arefiev 2005). 

However, there are observational evidences of the large 
thickness of the disk in SS433. The approximate equality 
of the amplitudes of the precessional modulation and of 
the primary eclipses (the latter measured in the preces- 
sional phases when the disk is maximum open) in opti- 
cal bands UFabrika fc Irsmambetova 20021 IFabrika 2004 |l 
indicates that the size of the donor is almost equal to 
that of the geometrically thick disk. The same relationship 
was found in different X-rays bands (Cherepashchuk et al. 
2004): from 1.5 keV to 100 keV both the precessional and 
eclipsing amplitudes increase while keeping comparable to 
each other. This means that for all acceptable binary mass 
ratios in SS433 {q = M^/Mop = 0.2 - 0.5) both the donor 
and the outer disk rim sizes are about ~ 0.4 — 0.5 in units 
of the binary separation and the estimate H/R ^ 1 looks 



Note here that the ratio of orbital to the break frequency of 
the order of 5-10 may indicate that the mass ratio in the binary 
system SS433 is larger than q > 0.3 (see results of Gilfanov & 
Arefiev 2005) 
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quite plausible. This may be a point of concern. For a disk 
with H/R ~ 1 the long viscous time ~ lOTV (recall that 
Tk relates to the Keplerian time at the outer disk radius) 
seems to be unlikely unless a small effective viscosity is 
assumed. 

So to have a long viscous time scale in the accretion 
disk in SS433 wc can admit that cither some outer parts 
of the disk (through which matter diffuses on the longest 
time scales) have a moderate thickness H/R ~ 0.2 — 0.3, 
or the viscosity parameter a in the accretion disk of SS433 
is smaller than 0.1 

The viscous time scale in an accretion disk can be 
treated as the time it takes for the accreting material 
to pass through the disk. This time may be estimated 
from the analysis of the nutational variability of jets in 
SS433 (see e.g . lKatz et al. 1 9821 ITHiTTis Newsom 19861 
IFabrika 2004 |l . The idea is that the nodding motions of 
the jets and the inner accretion disk have to be delayed 
with respect to tidal perturbations of the outer parts of 
the disk induced by the secondary companion. This time- 
delay was estimated to be of the order of a day. This means 
that the time for the passage of material through the disk 
is either extremely short (practically the free-fall time) or 
it is an integer multiple of the nutational period because 
of the periodic nature of the nodding motions. 

A short travel time of the material passing through 
the disk is in strong disagreement with our estimates in- 
ferred from the noise power spectrum. Wc should note 
here that the estimate of the travel time from the nuta- 
tional variability actually measures the propagation time 
of tidal perturbations from the outer parts of the accre- 
tion disk to the inner parts. Periodic tidal forces from the 
optical component may form a standing (in the binary's 
reference frame) structures in the outer accretion disk (see 
e.g. model calculations in Blondin 2000land references to 
earlier works therein) which can translate the nutational 
perturbations from outer parts of the disk to its inner 
parts on time scales much shorter than the matter travel 
time through the disk. In this picture, tidal perturbations 
may propagate through the disk with velocities close to 
the sound velocity. 

Of course, the disk in SS433 is far from being the 
standard one. It is strongly supercritical and is inclined 
to the orbital plane (w 20°). However, we can use as a 
guide the analogy with the twisted tilted subcritical ac- 
cretion disk observed in accreting binary system Her X-1. 
In that system the observed X-ray flux is strongly modu- 
lated with a period of ~ 35 days. The binary inclination in 
Her X-1 is close to 90 degrees, so we observe the system 
nearly edge-on. Over the 35-day period the X-ray emis- 
sion appears twice - the bright "main-on" state and then 
a weaker "low-on" state. In Her X-1 the observed 35-day 
modulation can be understood as the precessional motion 
of the outer parts of the inclined accretion disk (Gerend 
& Boynton 1976). The analysis of the beginning of an X- 
ray cycle (the so-called turn-on phase) observed by RXTE 
(Kuster et al. 2005) confirms that it is the outer edge of 
the accretion disk that opens the central X-ray source at 



the turn-on. The disk is not flat, it is warped (or twisted) 
because of the strong interaction of the innermost parts 
with the rotating magnetosphere of the central neutron 
star in Her X-1. The inner parts of the disk start screen- 
ing the central X-ray source at the end of the main-on 
X-ray state (see e.g. Shakura et al. 1999 and references 
therein). The disk in Her X-1 precesscs as a whole in the 
direction opposite to the orbital revolution, so the inner 
parts of the disk appear to lead the outer parts in the pre- 
cessional motion. This can be realized only when the tidal 
interaction from the outer parts is translated to the inner 
parts of the disk on a time scale which is shorter than the 
viscous time of the disk (of order of 20-30 days). 

3.3. Accretion disk imprints in SS433 variability 

Optical. It can be considered well established observation- 
ally that most of the optical variability in SS433 comes 
from the so-called accretion disk funnel, the open cone in 
the innermost part of a supercritical accretion disk (see 
e.g. IFabrika 200'4l for a review and recent calculations by 
Okuda et al. 2004). The funnel itself may by partially hid- 
den from the direct view by the strong wind from inner 
accretion disk. In this case some fraction of the optical 
radiation and variability will come from the wind photo- 
sphere. Both the funnel and the inner wind surround the 
bases of the jets. The structure of this region is not a sub- 
ject of this paper, the most important thing here is that 
the energy emitted in optical and UV bands is very large, 
of the order of ~ 10'**' erg/s. This undoubtedly tells us 
that most of the optical emission should be originating 
from reprocession of the internal X-rays generated in the 
inner accretion disk. Therefore in the framework of self- 
similar variability of the mass accretion rate in the disk 
( |Lyubarskii 1997| IChurazov et al. 2001(1 the optical vari- 
ability of SS433 should share the same time properties as 
observed in X-ray range. It naturally should have the self- 
similar power spectrum from / ^ t^j^^ to tj^^, where fvisc 
is the viscous time at the outer boundary of the accretion 
disk (the longest time scale in the disk) and tic is the light 
crossing time of the funnel. At frequencies above / ^ tj"^ 
the optical photometric variability should be smeared out. 

In SS433 the size (height) of the funnel can be evalu- 
ated to be of the order of 10*^ cm l|van den Heuvel 19811 
Marshall, Canizares, fc Schulz 2002| 

Revnivtse^'et'^L'^004ji and therefore the limiting 
frequency / ~ tj~^ ~ 0.01 Hz. Unfortunately, the existing 
data do not allow us to detect the decline of the power 
spectrum at frequencies higher than ~ 0.01 Hz (see Fig. 
^ but future observations with higher statistics may 
reveal it. 

It is interesting to note here that the variability of 
optical emission lines, which are emitted in the jets at 
distances above ~ lO*** cm (Borisov & Fabrika 1987) from 
the compact object where the jet cross size is of the same 
order, should start washing out at the jet light-crossing 
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frequencies of the order of ~ 10 ^ Hz, which should be 
detectable. 

X-ray. The X-ray flux of SS433 originates at the base 
of hot outflowing jet formed in the innermost regions of 
the accretion flow. Therefore we can naturally anticipate 
the X-ray flux to contain information about the mass ac- 
cretion rate variations present in the innermost parts of 
the disk which arc presumably brought in from the outer 
parts of the disk with the accreting matter. In this picture 
the similarity of the power spectra of variability of SS433 
in optical and X-ray spectral band is quite natural. Indeed, 
the direct correlation of variability of SS433 in optical and 
X-ray bands was recently discovered by Revnivtsev et al. 
(2004). 

Note that the X-ray variability should also be smeared 
out at high frequencies (in X-rays we observe the cool- 
ing jets). The limiting frequencies depend directly on the 
X-ray jets length r^, which presumably is ~ 10^^ cm 
IjWatson et al. 1986L Kotani et al. 1996). So a decline 
in the X-ray power spectrum at frequencies higher than 
0.01 — 0.1 Hz is expected. 
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Fig. 2. The cross-correlation of radio ('^2 GHz) and op- 
tical (photometric) variability. The light curves were re- 
binned into 2 X 10^ sec time bins. The positive lag means 
that the radio flux is delayed with respect to the opti- 
cal one. The shaded region denotes the rms uncertainties 
in the measurements of the cross-correlation at any given 
time lag (estimated from the data). The solid line shows 
the average cross-correlation function. The dashed line 
represents the average cross-correlation measured from the 
light curves rebinned into 10^ sec bins 

Radio. Most radio emission observed from SS433 orig- 
inates in the jets. The power spectrum obtained by us 
at very low Fourier frequencies (< 10~^ Hz, see Fig. 
clearly indicates that variability of radio emission in SS433 
has the same origin as the optical and X-ray variability. 
Even more direct confirmation of this hypothesis can be 
extracted by cross-correlating the optical and radio time 
series. From Fig. 12 we can see that the light curves at 
these spectral bands are significantly correlated. The ra- 
dio flux is delayed with respect to the optical photometric 



variability with a time lag of '^2 days. This time lag ap- 
proximately corresponds to the jet travel time from the 
accretion disk funnel (where most photometric variability 
is generated, see above) to the distance where the radio 
emission appears, ^ 10^^ cm (Paragi et al. 1999). 

The observed power spectrum (Fig. 1) indicates that 
there is no additional generation of the radio flux variabil- 
ity in jets on longer time scales (> 10^ sec). The points on 
the power spectrum at frequencies above ^ 10^^ Hz show 
that radio flux have a smaller variability amplitude than 
that found in the optical and X-rays at same frequencies. 
This can be a result of strong smearing of any possible in- 
trinsic variability in the radio flux on time scales shorter 
than the light-crossing time of the radio emitting region. 

4. Discussion 

As we have shown above, the straightforward interpreta- 
tion of the obtained broad-band power spectra of aperi- 
odic flux variations in SS433 can be made in terms of self- 
similar accretion rate modulations by the disk. How solid 
is this interpretation? It is based on the assumption that 
the observed frequency of the flattening of the power spec- 
trum can be associated with the viscous time scale of the 
accretion disk. This assumption seems to be supported by 
the analysis of power spectra of a dozen of LMXBs made 
by Gilfanov and Arefiev (2005). 

The accretion disk in SS433 is a complex structure 
and clearly is different from disks in LMXBs. A time- 
dependent hydrodynamic model of an accretion disk in 
a close binary (Blondin 2000) shows that the tidal forces 
truncate the accretion disk to radii of order half the aver- 
age radius of the Roche lobe and an effective a-parameter 
is of order 0.1 near the outer edge of the disk. The disk in 
SS433 has to be even more complex, it is inclined by an 
angle of 20 degrees to the orbital plane. So the resulting 
thickness of the disk must be further increased in compar- 
ison with any estimates based on the standard accretion 
disk theory. 

From the observational point of view there are clear 
indications that the inner disk around compact object in 
SS433 should be very thick {H/R ^1). If the a parameter 
of viscosity in such a disk is a ~ 0.1 the disk viscous time 
scale can not be as large as ~100 days. To obtain such a 
long viscous time, which are necessary to ensure the long 
time scale modulations of the mass accretion rate, we can 
assume that there is some (probably small) part of the 
outer disk with a moderate disk thickness H/R ^ 0.1 — 0.2 

A possible alternative explanation can be delineated as 
follows. Let the viscous time of the disk be actually small, 
of order of Tr- a few days. The short viscous time scale is 
favored from the nutational periodicity analysis. The disk 
will generate the power-law modulations in the accretion 
rate all the way down to the frequency 1/Tr- ^-^ 10^^ — 10^^ 
Hz irrespective of what kind of fluctuations are originally 
supplied with the incoming matter from the companion. 
So to explain the same power-law continuation of the spec- 
trum to lower frequencies, as observed, we need to invoke 
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some physical mechanism of generation of the power-law 
fluctuation spectrum from the donor star. In principle, 
such a mechanism could be a feed-back between the op- 
tical/UV flux fluctuations in the luminous central source 
(and the disk wind) and the illuminated photosphere of 
the donor, which will be transferred to the outer disk re- 
gion with matter passing through the vicinity of the in- 
ner Lagrangian point in the binary system. Indications of 
a strong heating effect in the donor of SS433 have been 
recently found through optical spectroscopy of the com- 
panion A-supcrgiant carried out at the processional phases 
close to the maximum disk opening (Cherepashchuk et al. 
2004). Depending on orbital and processional orientations 
(phases) different parts of the donor's surface should be 
exposed for the heating. 

The self-similar behavior of fluctuations is quite com- 
mon in various natural phenomena, and the power-law 
noise spectrum can be generated in the turbulent disk out- 
flow (e.g. Zeleny & Milovanov 2004 for a recent review). In 
this picture, the lower limit frequency of correlated fluctu- 
ations can be associated with star's precession frequency 
fpr ~ l/160(i ~ 6 X 10"® Hz, which does not contradict to 
the observed flattening frequency. However, the inspection 
of the obtained power spectrum does not reveal a signifi- 
cant break at frequencies about 10"'"' Hz corresponding to 
the short viscous time scale of the disk. 

So the simplest approximation of the obtained fluctua- 
tion power spectrum as one power law from 10"'' to 10"^ 
Hz appears to support the long viscous time scale of the 
accretion disk in SS433 and the interpretation of this spec- 
trum in terms of self-similar modulations of the accretion 
rate in the disk. 

5. Summary 

We constructed the broad band power spectra of aperi- 
odic variability of high-mass X-ray binary system SS433 
in radio, optical and X-ray spectral bands. The obtained 
properties of the spectra can be summarized as follows: 

— Power spectra of variability of SS433 in radio, optical 
and X-ray spectral bands have the power law shape 
(P DC /-") from ~ 0.1 Hz down to 10"'^ Hz with flat- 
tening afterwards. At frequencies > 10"^ Hz the power 
law index is approximately a ~ 1.5 

— The power spectrum of radio variability presumably 
has one more break (steepening) at frequencies ^ 
10~^ Hz which can be a result of smearing with large 
light-crossing time of the radio emitting region 10^ 
sec) 

— Radio and optical (photometric) fluxes of SS433 are 
correlated and the radio flux is delayed by about 2 
days with respect to the optical one. The correlated 
variability in optical and X-ray spectral bands was pre- 
viously detected by Revnivtsev et al. (2004) 

The observed behavior of fluxes of SS433 in different 
spectral bands can be interpreted in the framework of self- 
similar accretion disk variations proposed by Lyubarskii 



(1997) and elaborated by Churazov et al. (2001) and 
Gilfanov & Arefiev (2005). The accretion disk introduces 
the mass accretion rate variations at any given radius for 
example due to fluctuations of the viscosity parameter a. 
The resulted variable mass accretion rate enters the re- 
gion of the maximum energy release near the compact 
object and generates powerful variations in the emitted 
X-rays which then are reradiated in optical and UV spec- 
tral bands by a funnel surrounding the jets. We do not see 
directly this X-ray emission from the innermost parts of 
the accretion flow because of the high inclination of the 
binary system and a large optical depth of the geometri- 
cally thick accretion disk and the outflowing wind. Mass 
accretion rate fluctuations are translated to the density 
variations in the outflowing hot baryonic jet and give rise 
to the observed X-ray variability correlating with the op- 
tical one (Revnivtsev et al. 2004). As the expanding jet 
cools down, it starts emitting in optical lines and in ra- 
dio. The time lag between the radio flux and the optical 
photometric flux (emission of the accretion disk funnel) 
equals to the jet travel time from the funnel to the radio 
emission region. 
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